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We present the first constraints on the full parameter space of the Galileon modified gravity model, consider- 
ing both the cosmological parameters and the coefficients which specify the additional terms in the Lagrangian 
due to the Galileon field, which we call the Galileon parameters. We use the latest cosmic microwave back- 
ground measurements, along with distance measurements from supemovae and baryonic acoustic oscillations, 
performing a Monte Carlo Markov Chain exploration of the 9-dimensional parameter space. The integrated 
Sachs-Wolfe signal can be very different in Galileon models compared to standard gravity, making it essential 
to use the full CMB data rather than the CMB distance priors. We demonstrate that meaningful constraints are 
only possible in the Galileon parameter space after taking advantage of a scaling degeneracy. We find that the 
Galileon model can fit the WMAP 9-year results better than the standard A-Cold Dark Matter model, but gives a 
slightly worse fit overall once lower redshift distance measurements are included. The best-fitting cosmological 
parameters (e.g. matter density, scalar spectral index, fluctuation amplitude) can differ by more than 2cr in the 
Galileon model compared with ACDM. We highlight other potential constraints of the Galileon model using 
galaxy clustering and weak lensing measurements. 



I. INTRODUCTION 

One of the major unsolved mysteries in cosmology is 
the nature of the 'dark energy' that is causing the observed 
present-day accelerated expansion of the Universe lil nlOJ . In 
the context of General Relativity (GR), the vacuum energy, 
also known as the cosmological constant. A, is the simplest 
explanation for dark energy. In the current standard cosmol- 
ogy, known as the A - Cold Dark Matter (ACDM) model, 
the cosmological constant accounts for approximately 70% 
of the total energy content of the Universe today, with non- 
relativistic dark and baryonic matter, neutrinos and photons 
completing the cosmic inventory. ACDM performs very well 
against observations. However, this model has serious theoret- 
ical problems, the biggest one being the huge discrepancy of 
many orders of magnitude between the standard quantum field 
prediction for A and the value inferred from observations. It 
has become widely accepted that one needs to look beyond A. 
A possible and attractive alternative is to modify GR on large 
scales in such a way to allow cosmic acceleration without the 
need to invoke a cosmological constant or any other negative 
pressure fluid iTTl [T2l . Interest in modified gravity models 
has been growing over the past decade [131, with significant 
progress being made in both analytical modelling L14h22< and 
numerical simulations |23 -25|. 

Long before the discovery of the accelerated expansion, in 
1974, Horndeski derived the most general action for a single 
scalar field that yields equations of motion (EOM) which con- 
tain up to second order derivatives of the fields [26 1, avoiding, 
therefore, the presence of Ostrogradsky ghost degrees of free- 
dom ll27l . However, it was only recently that such an action 
started to gain relevance in cosmology Il28ll29]| . A large num- 
ber of successful modified gravity models can be considered 
as special cases of the general Horndeski theory, which satisfy 
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certain desired properties or symmetries. The most famous 
recent examples include the Galileon model 1.30.- 321, kinetic 
gravity braiding I33ti35l . Fab-Four ||36H39]| . Fab-Five HO), k- 
mouflage |41| and others P2ti45]| . 

In this paper we focus on the Galileon model, which was 
originally derived by generalizing the four-dimensional effec- 
tive boundary action of the braneworld DGP model B6H491 . 
In both these models, the modifications to gravity are deter- 
mined by an extra scalar degree of freedom if, known as the 
Galileon, whose Lagrangian is invariant under the Galilean 
shift symmetry d^ip — ^ d^ip + b^, where is a constant 
vector. The branch of the DGP model that leads to cos- 
mic acceleration suffers from ghosts ||50H52]| : the Galileon is 
free from such problems. In a four-dimensional Minkowski 
space there are only five Galilean invariant Lagrangians that 
are second order in the equations of motion, despite contain- 
ing highly nonlinear derivative couplings of the scalar field 
1301. In order to investigate the cosmological implications of 
the Galileon model, its Lagrangian needs to be generalized 
to curved spacetimes and in 1311 132| the authors have shown 
that explicit couplings between the Galileon field derivatives 
and the curvature tensors are needed to retain the equations of 
motion up to second order, although they explicitly break the 
Galilean symmetry. 

The presence of the nonlinear derivative interactions and 
the couplings to the curvature tensors in the Galileon model 
change the way in which matter determines the geometry of 
spacetime relative to the GR prediction. However, near mas- 
sive objects like the Sun, the laws of gravity are constrained to 
be very close to GR and, therefore, any viable modified grav- 
ity model needs to revert to the 'standard' solution on these 
scales f5S'-'55l. Remarkably, the nonlinear derivative terms al- 
low, at the same time, for the suppression (or screening) of the 
deviations from GR by effectively decoupling the scalar field 
from gravity in high curvature regions by a mechanism known 
as Vainshtein screening |56 1. In brief, near matter sources, the 
nonlinear terms of the equations of motion become important, 
strongly suppressing the spatial gradient of the scalar field. 
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which is the extra 'fifth force', on scales below the so-called 
'Vainshtein radius'. On the other hand, beyond the Vainshtein 
radius the nonlinear terms are subdominant, such that both 
the Galileon field and the Newtonian potential satisfy a linear 
Poisson equation, and the total gravitational force gets a non- 
negligible contribution from the 'fifth' and the 'normal' GR 
forces. The Vainshtein mechanism is crucial to the success of 
the Galileon model and it is in many respects similar to the 
implementation of other screening mechanisms such as the 
chameleon El El, dilaton fHllllEll, symmetron Il60ti63l 
and disformal screening |64|. The main difference is that in 
the case of the Vainshtein mechanism, the range of the screen- 
ing depends only on the properties of the gravitational source 
as opposed to the chameleon mechanism, in which case the 
screening depends also on the cosmological environment. 

The way the background cosmological evolution can be af- 
fected in Galileon models has akeady been studied and con- 
strained using observational probes of the geometry of the 
Universe such as type la Supernovae (SNIa), the scale of the 
Baryonic Acoustic Oscillations (BAO) feature in the galaxy 
distribution and the position of the peaks of the Cosmic Mi- 
crowave Background (CMB) angular power spectrum of tem- 
perature fluctuations ll65lj7Tl . Furthermore, by studying the 
evolution of density fluctuations in linear perturbation the- 
ory, it has also been shown how Galileon models can be con- 
strained theoretically, by the imposition of conditions to avoid 
the presence of ghosts and other pathologies 1671 |69l , and ob- 
servationally, by investigating the predictions for the growth 
rate of Large Scale Structure (LSS) |70--73'] and for the fuU 
CMB, linear matter and weak lensing potential power spectra 
ll74l (see also ll75l for a study of the bispectrum of the matter 
fluctuations). These studies revealed a number of distinctive 
features of the Galileon model, which could, in principle, be 
used to place strong constraints on the model and distinguish 
it from the ACDM paradigm and from other modified gravity 
theories. 

In Il74l we took the first step towards the thorough explo- 
ration of the Galileon parameter space that we present in this 
paper In particular, by using a modified version of the pub- 
licly available CAMB code |76|, we showed that the Galileon 
model predictions depend sensitively on the values chosen for 
the model parameters. We demonstrated that the Integrated 
Sachs-Wolfe (ISW) effect can play a key role in constraining 
the Galileon parameter space, thanks to the distinctive signa- 
ture the model has on the largest angular scales of the CMB 
power spectrum. Our previous work has therefore motivated 
the analysis presented in this paper, where we explore the full 
cosmological parameter space and not only its Galilean sub- 
space, using Monte Carlo Markov Chain methods (MCMC) 
with the aid of the publicly available CosmoMC code [77|. To 
allow the rest of the cosmological parameters to vary is crucial 
to fully understand the degeneracies that might exist between 
them, and it is, ultimately, the only way to truly understand 
how well a given model is able to reproduce the observational 
data. 

The present paper is structured as follows. In Section[lI|we 
introduce the covariant Galileon model, present the fully co- 
variant and gauge invariant perturbed field equations we de- 



rived in f74] and that we use to obtain the results. We also 
recap the main physical features of the model. In Section 
[HI] we specify and discuss the full cosmological parameter 
space which we constrain, the theoretical priors that determine 
the viability of some regions of the parameter space and the 
datasets that we use to derive the constraints. In particular, we 
point out that if all of the parameters that enter the model La- 
grangian are allowed to vary, then there is an infinite degen- 
eracy region along the Galilean subset of the full parameter 
space. We present and analyse the results in Section[rV| where 
we look at how different combinations of data constrain the 
Galileon model, highlighting the tensions between the best- 
fitting models derived from different datasets. We present our 
conclusions in Section IV] 

Throughout this paper we assume the metric convention 
(+,—.—,—) and work in units in which the speed of light 
c = 1. Greek indices run over 0, 1, 2, 3 and we use SttG — 
K — Mp^ interchangeably, where G is Newton's constant and 
Mpi is the reduced Planck mass. 



11. THE GALILEON MODEL 

In this section we briefly review the covariant Galileon 
model and the field equations used to obtain the physical pre- 
dictions of the model. For a detailed derivation we refer the 
reader to |74| (see also [[78l ). 

The full action of the covariant Galileon model, which has 
no direct coupling between matter and the Galileon field (see 
however [69..74.,79-.--84J ). is given by 
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where g is the determinant of the metric g^^, R is the Ricci 
scalar and the model parameters ci_5 are dimensionless con- 
stants. The five terms in the Lagrangian density, fixed by the 
Galilean invariance in flat spacetime, d^ip d^Lp + h^, are 
given by 



c. = 



Vp</jV^</j 2(0(^)2 -2(V^V,(^)(V^VV) 



1 

-RV^ipV''ip/2 

1 

M9 



- 3(nvp)(v^v,(^)(V'^v» 



+2(V^VV)(V.V''(^)(VpV''(p) 
-6(V^(^)(V^VV)(V''(^)G,p 



(2) 



in which (p is the Galileon scalar field and = Mp\Hq 
with Hq being the present-day Hubble expansion rate. The 
derivative couplings to the Ricci scalar R and the Einstein 
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tensor G^,y in £4 and £5 are necessary to prevent the EOMs 
from having higher than second-order derivatives of the metric 
and the Galileon field in curved spacetimes, such as the one 
described by the Friedman-Robertson- Walker (FRW) metric 
f31J . Such terms, however, break the Galilean shift symme- 
try. 

As is customary in Galileon studies, one is generally inter- 
ested in cases where cosmic acceleration is due only to the 
field kinetic terms and therefore we will set the potential term 
ci to zero. Note also that an important difference between the 
Galileon model and other subsets of Horndeski theory is that 
there are no free fuctions in the action of Eq. ([T]). In particular, 
this prevents the Galileon model from having a ACDM limit 
which in the end plays an important role in distinguishing be- 
tween these two models (see however Il85ll86l ). 

The modified Einstein equations and the Galileon EOM are 
obtained by varying the action of Eq. ([T]), with respect to g^i, 
and (y3, respectively. Given the length of these equations we do 
not show them in this paper. Interested readers can find them 
in Eqs. (A1-A7) of Ref. 1741 . 



A. Friedmann equations 

In the Galileon model, the structure of the two Friedmann 
equations is not changed with respect to standard GR. In this 
paper we will always consider the case of a spatially flat Uni- 
verse for which the first and second Friedmann equations are 
given by: 



where and £7^ are the present-day values of the fractional 
energy density of the radiation and matter components, re- 
spectively. The background pressure of the radiation and mat- 
ter is given, respectively by. 
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B. Galileon field equations 

To describe the Galileon field up to linear order in pertur- 
bations, we shall make 3 + 1 spacetime decompositions of the 
tensor quantities ll87l [88l . This is achieved by using the pro- 
jection tensor h^^j, which is defined as h^^ = g^^ — u^Ui, and 
can be used to obtain projected covariant tensors which live in 
3-dimensional hyperspaces perpendicular to an observer's 4- 
velocity u^. For example, the covariant spatial derivative V 
of a tensor field T' 
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is defined as 
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The energy-momentum tensor and the covariant derivative 
of the observer's 4-velocity field can be decomposed, re- 
spectively, as 



Tf^u ^ T^tiu + 2q(^f_iU^) + pUf_,u^ - ph^^, (8) 

,9, 
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where 6 = 3a/ a = 3H, where a is the cosmological scale 
factor, H the Hubble expansion rate (throughout the paper we 
shall use 6 and H interchangeably); the overbar denotes back- 
ground averaged quantities, the overdot denotes a derivative 
with respect to cosmic time and we are preserving the nota- 
tion used in 1741 1 . 

In Eqs. ^ and Q, p = +Pm +Pip and pr + pm + p^ 
are, respectively, the total pressure and energy density in the 
Universe. The subscripts r, m and ^ denote the contribution 
from the relativistic degrees of freedom (photons and mass- 
less neutrinos), non-relativistic matter (dark and baryonic) and 
the Galileon field, respectively. In the case of the uncoupled 
Galileon model, the cosmic evolution of pr and is fixed by 
their present-day values p^o and p,no, as 



pr = Pro a 
Pm — PmOO' 



where tt^^ is the projected symmetric and trace-free (PSTF) 
anisotropic stress, is the heat flux vector, p is the isotropic 
pressure, p is the energy density, ct^^ the PSTF shear tensor, 
'^f.iu — the vorticity, 6 = V"Wq = 3H as intro- 

duced above and — ii^ is the observer's 4-acceleration. 
The time derivative can be expressed in terms of the covariant 
derivatives as = u"Vq,0. Square brackets mean antisym- 
metrization and parentheses symmetrization. The normaliza- 
tion is u"ua = 1, in accordance with our choice for the metric 
signature. 

The physical quantities in the total energy-momentum ten- 
sor T^y, Eq. (|8]l, have contributions from each of the matter 
species in the universe: 



p = p +p 

P = P'' + P'" 
9^ = 9m + 9^ 
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In ll74l we presented a detailed explanation of how to derive 
the Galileon field contribution to the quantities of Eqs. ( [T0| - 
13 I. However, for brevity, in this paper we limit ourselves to 
quoting the final result, which, up to first order in perturbed 
quantities, is given by 
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in which □ = V'^V^. In Eq. ([itJ, ^ l^^SJ v)'^ and 
£iiv = '^°'u^y^ iiav(i are PSTF rank-2 tensors, both of which 
live in the 3-dimensional hypersurface perpendicular to 
(u^V t^i^^S/^-^ip = u^S/^i, = 0), where W^q^^ is the Weyl cur- 
vature tensor. The equality symbol = means that we have ne- 
glected the terms which are higher than linear order in small 
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perturbations. Angular brackets indicate trace-free quantities. 

Note that the £2 and £3 terms do not contribute to the 
Galileon field anisotropic stress. 

The Galileon field EOM, again, up to linear order in small 
perturbations, is given by 
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Eqs. (14 - 18 I also contain the background (zeroth-order 



terms) as the non-hatted terms. Here, for clarity of exposition 
and because these terms enter the physical calculations sepa- 
rately, we extract the background expression for the Galileon 
field density, pressure and equation of motion, which follow, 
respectively, as: 
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where we have neglected the overbars on Lp to lighten the no- 
tation. The heat flux q'^ and anisotropic stress tt^^ contain 
only perturbed quantities and vanish at the cosmological back- 
ground level. Note that assuming the de Sitter limit where 6 
and ip are both constants, it can be confirmed that Eqs. ( 19 1 
and ([20| satisfy p^p + = 0, that is w = —1. 



Eqs. g |4] |5] |6] [19] |2U] |2T) form a set of differential equa- 
tions for the background evolution of the Universe, whereas 
Eqs. ([T4|-[T8]l are the expressions which enter the equations 
in CAMB to solve for the evolution of the linear perturbations 
01761. 



C. Cosmology of the Galileon model 

The cosmological features of the Galileon field have been 
studied in the literature in great detail both at the background 
116514671 170 78] and at the Unear perturbation level |69 71 74 
l89l . Here we just briefly outline the most distinctive aspects 
of the covariant Galileon model. 



the quasi-static limit (the limit where the time derivatives are 
neglected relative to spatial derivatives) is in general a good 
approximation for the full equations. In [91] the authors de- 
rived the linear perturbation equations for the full Horndeski 
Lagrangian, which includes the Galileon model studied here. 

In |[74l we showed that the Galileon field can cluster 
strongly at the linear level since early times and in a way 
comparable to the clustering of matter Such a strong cluster- 
ing signal can be attributed to the non-trivial behavior of the 
pressure perturbation and the anisotropic stress in the Galileon 
model. 

The clustering of the Galileon field is intimately related 
with the time evolution of the gravitational potential. In the 
Galileon model, the gravitational potential can evolve signifi- 
cantly even deep inside the matter dominated era (in contrast 
to the ACDM paradigm, where it remains constant throughout 
matter domination). This leads to very clear signatures in the 
low-Z region of the CMB temperature power spectrum where 
the ISW effect, precisely determined by the time evolution of 
the gravitational potential, makes the dominant contribution. 
In addition to a nonstandard time variation, the gravitational 
potential can experience an overall deepening with respect to 
the ACDM prediction, which also leaves distinct imprints on 
certain cosmological observables. In particular, it can strongly 
modify the growth rate of matter density fluctuations ll7TH74ll 
and also make very different predictions for observables such 
as weak lensing and the cross correlation between the ISW 
effect and the galaxy distribution 11741 . 



III. METHODOLOGY 



1. Background 

In the authors have shown the existence of a tracker 
solution of the dynamical background equations which ap- 
proaches de Sitter evolution {w — —1) at the present time. 
This happens after periods of radiation and matter domination, 
thus allowing for a viable expansion history. The tracker is 
characterized by phantom evolution (w < —1) of the Galileon 
field in the past and works as an attractor for other solutions 
with different initial conditions. In iTTSlI it was shown that 
the background cosmological data prefers solutions which ap- 
proach the tracker at late times so that the phantom period 
of the evolution only occurs close to today and the equation 
of state parameter w does not become too much smaller than 
w = —1. 

A noteworthy aspect here is the possibility of having ghost- 
free phantom dynamics ll 22l [33l [34l which could lead to clear 
signatures in the expansion rate of the Universe I.90J , thus 
helping to distinguish this model from other models such as 
ACDM or Quintessence. 

2. Linear perturbations 

The full linear perturbation equations of the Galileon model 
have been derived in 11741 |89l , where it was also shown that 



A. The Galileon cosmological parameter space 

We restrict ourselves to a flat geometry for the background 
space-time. In this case the full cosmological parameter space 
we consider is eleven-dimensional with six cosmological pa- 
rameters: f2°, r2°, h, T, Us and log[10^''As], which are, re- 
spectively, the present-day values of the fractional densities 
of dark matter and baryonic matter, the Hubble expansion rate 
h — iJo/(100 km/s/Mpc), the optical depth to reionization, 
the scalar spectral index (n^) and amplitude (As) of the pri- 
mordial power spectrum of the scalar fluctuations; plus the 
five Galilean parameters C2, C3, C4, C5 and (fi, where (pi is 
background value of the time derivative of the Galileon field at 
the initial time ti (orredshift Zi) when the calculation starts. In 
this paper we take — 10^. The energy density of radiation is 
held fixed at ilro^i^ = 4.18 x 10^^, which is given by the tem- 
perature of the CMB photons, Tcmb = 2.725K, with an ef- 
fective number of massless neutrinos given by ]\[™^^^<^^^ ^ — 
3.04 |7|. 

Since we are assuming a spatially flat background Universe, 
we require that the fractional energy densities today satisfy 
fi" + + + = 1- This constraint equation means 
that one of the Galileon parameters C2-5, (pi can be uniquely 
derived from the others. In this work, we choose C2 as the de- 
rived parameter, and follow a trial-and-error approach to find 
the value of C2 for which the universe is spatially flat. Note 
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that for each new trial value of C2, we need to evolve the whole 
set of background equations from z = to ^ = 0, to check 
if the constraint is satisfied. This procedure can take a signifi- 
cant time, particularly because we have to sample many points 
in the parameter space during the MCMC search. To speed 
up this process, we designed an algorithm that takes adaptive 
step lengths in setting the trial value of C2, as well as adopting 
a conservative criterion by stopping the trial-and-error search 
when \n"^ + n"^ + n° + fl°-l\ < 10-^. This accuracy is 
sufficient for all our numerical results. Because C2 is now a 
derived parameter, the dimensionality of the parameter space 
is reduced by one. 

The value of the Galileon background density j5tp_i, at the 
starting redshift, is determined by the values of (fi and 9i, 
the latter being given by the fixed matter and radiation com- 
ponents via Eqs. (|3] |5]l (the Galileon energy density is neg- 
ligible at the very early times). Because this is a more 
physically meaningful parameter than (pi, despite sampling 
through different values of (pi in the MCMC search we use 
the ratio between the Galileon and matter energy densities 
at the starting redshift, Pip.i/pm,i, when we quote the ob- 
servational constraints below. Moreover, Eq. ( [T9| can be 
used to translate in between the values of p^^i and (pi via 
Pip,i ~ '''csi^f 6'?/(9Af^). This approximation holds at suf- 
ficiently early times where the Galileon density is negligibl^ 
provided also that C5 is not too much smaller than the values 
of the other c„ (which is the case for the best-fitting regions 
of the parameter space). 

To solve for the evolution of the linear perturbations, we 
also need the initial conditions of the Galileon field perturba- 
tion and its time derivative at Zi, which are, in principle, also 
free parameters. However, we know that these numbers must 
be very small at high redshift. As a result, we have set both 
of them to be exactly zero, and checked that our results are 
insensitive to sufficiently small changes around these initial 
values for all length scales (or fc-modes) that are of interest to 
us. Here, 'sufficiently small' means small enough to still be in 
the regime of linear perturbation theory. Typically, such a con- 
dition is quoted as Sipi <C (pi- However, the background value 
of the field is irrelevant in the covariant Galileon model since 
the background equations only involve ip and (p and not (p. As 
an alternative, we adopt 6pi ^ (pi /Hi as a criterion for the 
validity of linear perturbation analysis, and this restricts Sipi 
and d(pi to be so small that it makes no practical difference if 
they are set to be exactly zero. 



B. Scaling degeneracy in tlie Galileon model 

By looking at Eqs. ([T4]-[T8|, which are the equations that 
fully govern the Galileon physics, we note that they are invari- 



ant under the following transformations: 
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= Ci/B^ 
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= C5/B^ 
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^ipB, 



(22) 

in which B is an arbitrary constant and the transformation of 
ip holds for both the background and perturbation parts. This 
scaling relation is what allows [67] and [7 1 ] to use different 
sets of parameters to characterize the Galileon model. The 
reason for this scaling lies in the fact that each of the Galileon 
Lagrangians £i yields terms which all have the same power 
in the Galileon field p (for the counting of the power, the time 
and spatial derivatives of ip are treated equally as ip). Note 



that the transformation of Eqs. (22 1 preserves the signs of the 
parameters C2 and C4, but not those of C3, C5 and (pi (if B < 
0). This illustrates that the scaling also contains a 'reflection 
symmetry' in the C3 — C5 — <pi subspace, which is associated 
with the sign of B, instead of its magnitude. 

The existence of such a scaling relation is non-trivial in 
dark energy and modified gravity models. For example, the 
EOM of a Quintessence field (p with a self-interacting po- 
tential V{(p) is ip + 9ip + dV{(p)/dip = 0; therefore, unless 
V{(p) cx ip'^, it is impossible to allow for this scaling relation. 
From a practical point of view, the scaling of the Galileon 
field (p is realized by rescaling its time derivative at the initial 
time, Pi^ As a result, according to Eqs. (22i, the impact of 



smaller values of (pi can always be compensated by larger val- 
ues of the c„ parameters and vice versa, thus making the c„ 
parameters unbounded and preventing a proper constraint of 
the parameter space. 

The scaling relation allows one to further reduce the dimen- 
sionality of the Galileon subspace of parameters by one. This 
can be done by using one of the Galileon parameters as a refer- 
ence to write down invariant quantities under the scaling (e.g., 
in 1 71 1, the reference parameter is the present-day value of the 
Galileon field time derivative). In this paper, we take C3 as the 
reference parameter and therefore the invariant quantities are 



C2 C4 



C5 1/3 



2/3' 4/3' 5/3' 



(23) 



Note that p,^.i/Pm,i is, by definition, invariant under the scal- 
ing (c.f. Eq. ([T4|). When running the Markov chains, while 
one could allow for all of the parameters to vary and apply 
the constraints to the set of Eq. ( [23] l, it is easier to fix the 
value of C3. This will increase the convergence rate of the 
MCMC algorithm since now there is one fewer dimension to 
sample from and the unboundedness of the c„ parameters can 
be avoided. 



' Otherwise, one would need to take the initial Galilean density into account 
in the determination of 8i, which is then given, approximately, by 6i 
y/3K{pmi + pri), via Eq. ja . 



^ In principle, to achieve the exact scaling at the linear perturbation level, one 
has to resize the Galileon field perturbation Sipi and its time derivative (5</3j, 
accordingly. However, thanks to our choices of initial conditions, namely 
Si^i = hip I = 0, such resizing does not need to be done explicitly. 
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Note that in principle any other Gahlean parameter could 
be used as the reference to write down the scaling invariant 
terms. In practice, however, C2 cannot be used as the fixed pa- 
rameter since it is the parameter which we tune to yield a con- 
sistent cosmological background evolution with the required 
amount of dark energy to make the Universe spatially flat, and 
is unknown a priori. Practically, it is not a good idea to fix C4 
or C5 either, because when all parameters are unfixed they are 
more likely to cross or become very close to zero - we have 
checked that, for example, fixing C5 = 1 (which is roughly 
the best-fitting value of C5 when C3 = 10) causes C3 and C4 
to become essentially unbounded since C5 can be as small as 
10^*^ when C3 — 10. However, having derived the constraints 
with C3 fixed, one can always use the scaling of Eqs. (|22]i to 
scale the constrained regions to the case where, for example, 
C2 = —1, which, in our convention in Eq. ([T]l, corresponds 
to the Galileon field having a standard scalar kinetic term, but 
with a different sign. 

In summary, the nine-dimensional parameter space we aim 
to constrain in this paper is specified by 

ncoMbQ,h,T,ns,\og [10^° A,] , -^^cl^^if^ , 
\ ^3 C3' / 

in which the first six are cosmological parameters and the 
last three are the Galileon (physical) parameters. The remain- 
ing Galileon parameters, namely 020^^^'^ and p,p^i/pm,,i, are 
treated as derived parameters. 

C. Theoretical constraints 

In addition to the observational constraints that we will out- 
line in the next section, the Galilean parameter space is a pri- 
ori constrained by the requirement to avoid the appearance 
of theoretical instabilities. The Galileon Lagrangian, being 
a subset of the more general Horndeski Lagrangian, is auto- 
matically protected against the propagation of Ostrogradsky 
ghosts as the equations are retained up to second order [27 1. 
However, other sorts of theoretical pathologies may still arise. 

In this paper we consider the conditions for each point in 
parameter space not to develop ghost degrees of freedom or 
Laplace instabilities in the scalar sector of the linear pertur- 
bations (see e.g. L69l [85l l92l for a discussion and derivation 
of the stability conditions). The no-ghost and no-Laplace sta- 
bility conditions, despite applying to the scalar perturbations, 
depend only on background quantities [69 1. When the MCMC 
algorithm tries a new point in parameter space, our code first 
solves the background evolution, testing whether or not it sat- 
isfies all the stability criteria. The calculation of the evolution 
of the perturbations and the subsequent likelihood evaluation 
is only performed if the point is theoretically viable. We only 
test the theoretical stability of any given point in the past, since 
there is no evidence that the instabiUties cannot develop in the 
unprobed future. 

One could also consider other theoretical conditions such as 
those which ensure that the Galileon field perturbation does 
not propagate superluminally, i.e. cj. > I (see e.g. ||93l ). 



However, such cases do not necessarily imply the existence 
of pathologies such as the violation of causality (see e.g. 
||94l [9511) and therefore we do not employ them. We also do 
not rule out a priori cases where < Oat some point in time, 
but instead let the data decide thek viability. 

We want to stress that the theoretical constraints are a con- 
venient way to select only those points which give viable per- 
turbation evolution, and once these constraints are satisfied so 
that a trial parameter point is not rejected straightaway, they 
play no further role in the calculation of likelihoods. 



D. Datasets 

To derive the constraints on the Galileon parameter space 
we use the CMB data of the full WMAP 9-year final release 
|7 1 in the form of the FORTRAN likelihood software provided 
by the WMAP team, for which we have adapted CosmoMC. 
Despite being WMAP's final release, the dataset will soon be 
replaced by the upcoming Planck satellite data. However, we 
expect that the differences between these two datasets will be 
of little importance as regards the Galileon model constraints 
because, as we will see in the following sections, the major 
impact of the Galileon model on the CMB happens on the 
largest angular scales, where Planck is not expected to per- 
form much better than WMAP due to the limits imposed by 
cosmic variance. Note that unlike ll68l 171117211 . we do not use 
the WMAP distance priors, which are derived parameters and 
require assumptions to be made about the background cosmol- 
ogy. Moreover, below we will see how the Galileon model can 
have an impact on the largest angular scales of the CMB and 
also on the amplitude of the acoustic peaks, which illustrates 
the advantage of using the full CMB data over using only the 
information encoded in the positions of the acoustic peaks. 

As a low redshift probe, we consider the three year sample 
of the Supernova Legacy Survey (SNLS) project [1| which 
contains 472 type la supemovae ranging from z « 0.15 to 
z w 1.1. When quoting the results, we marginalize over the 
parameters that are used in the calibration of the intrinsic lu- 
minosity of the SNIa of the SNLS sample. 

To complement the constraints from SNIa luminosity dis- 
tances, we also use measurements of the angular scale of the 
BAO feature in the galaxy distribution whose dependence on 
the background expansion rate differs from the SNIa con- 
straints and, therefore, probes a different region of the pa- 
rameter space. We use the BAO measurements from the 6df 
Galaxy Survey at redshift z = 0.106 |3 1, from the Sloan Dig- 
ital Sky Sui-vey (SDSS) DR7 at 2 = 0.275 and z = 0.35 
and from the SDSS-III Baryon Oscillation Spectroscopic Sur- 
vey (BOSS) at z = 0.57 0. 



1. Why we do not use growth rate and clustering data 

Several studies have shown that the modifications of grav- 
ity in the Galileon model can significantly enhance the growth 
of linear matter fluctuations on sub-horizon scales 1691 ITTl - 
|74l|89j. In particular, in t71n73J . the authors used measure- 
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ments of the linear growth factor to place strong constraints 
on the Galileon model. These studies have shown that there 
is an unavoidable tension in the Galileon model between its 
ability to fit, at the same time, background and growth data. 
However, these conclusions assume the validity of linear per- 
turbation theory in the Galileon model on the scales where 
the growth factor is measured. In GR, one would expect the 
scales probed by current growth factor measurements to be 
mildy nonlinear 1961 . However, for the Galileon model we 
do not yet have reliable knowledge about the true scale on 
which the Vainshtein screening starts to be important and how 
it affects structure formation: as this is a purely nonlinear ef- 
fect, it is by definition not present in a linear theory analy- 
sis. For example, numerical simulations have shown that in 
other modified gravity models such as the f{R) and dilaton 
159, 97, ,98] , linear perturbation theory can fail even on scales 
as large as A; ~ 0.01 /iMpc~^ because of the chameleon 
screening ll24l l25l l98l. It is reasonable to suspect that simi- 
lar situations might occur in the Galileon model. 

Our modifications to the CAMB code allow us to obtain the 
Galileon predictions for the linear matter power spectrum and 
growth rate ll74ll and we could in principle use these to place 
further constraints on the model. However, given the above 
reasoning we remain cautious about using clustering data for 
now. We argue that a better understanding of the true im- 
pact of the Vainshtein screening is needed before attempting a 
more rigorous confrontation of the predicted clustering power 
and growth rate with the observational data. Such an analy- 
sis of the Vainshtein effect in the Galileon model will be left 
for future work (see also |99 - 108 1 for some existing studies 
in this direction). In the present paper, we limit ourselves to 
what the CMB and background data can tell about the param- 
eter space of the Galileon gravity model. 



IV. RESULTS 
A. Numerical details 

We use the publicly available CosmoMC code ITTTI to carry 
out the formal exploration of the full Galileon cosmologi- 
cal parameter space, using the Metropolis-Hastings algorithm 
(see ITTTI for a concise description) as the MCMC method to 
draw the samples from the posterior probability distribution, 
which is determined by the observational data. The code is 
compatible with our modified version of the publicly available 
CAMB code [76 1 , which solves for the background and pertur- 
bation evolution in Galileon models |74 |. The results we show 
in this section were obtained by running six chains in parallel 
with the stoping convergence criterion R — 1 < 0.02, where 
R is the Gelman and Rubin statistic given by R ="variance of 
chains means'V'mean of chains variances" [109]. The MCMC 
algorithm is used to explore the parameter space with linear 
steps along all directions, except the (pi direction, where the 
steps are logaritmic since this parameter can take values that 
differ by a few orders of magnitude. The estimation of the 
likelihood from the samples was performed using the routines 
in the Getdist software supplied in the CosmoMC package. 



We do not consider the first half of the chains in the likelihood 
evaluation to eliminate points sampled during the 'burn-in' 
period of the chains. When plotting the marginalized parame- 
ter distributions and the respective probability limit contours, 
we smooth the likelihood surface by applying a Gaussian fil- 
ter (in addition to the smoothing already performed by the 
Getdist routines). We made sure that the smoothing does 
not affect our results. In Table |I] the results are for the un- 
smoothed distributions as determined by the Getdist eval- 
uation. We have also made sure that the prior range limits 
specified for each of the parameters were sufficiently far away 
from where the likelihood distribution is non-negligible. 

In addition to the standard modifications needed to link 
CAMB with CosmoMC when new parameters are added or 
alternative cosmological models are studied, we have also 
changed the way in which the h cosmological parameter is 
sampled. In the latest version of the code, the MCMC al- 
gorithm samples over the parameter rl""^ /(F^'^ (known in the 
code as the parameter thet a) where r^"^ and d^'^ are, respec- 
tively, the sound horizon at recombination and the angular di- 
ameter distance to recombination. Having sampled rl°'^ / d'-^'^ , 
the code then finds the value of by a trial-and-error proce- 
dure. The reason why it is done this way is because rl'^'^ / d'-^'^ 
is much less correlated than h with the rest of the cosmolog- 
ical parameters, thus improving the performance of the pa- 
rameter space exploration. However, in our case one has also 
the trial-and-error search for C2 which can be in conflict with 
the search for h. The value of C2 is fixed when the code is 
trying values of h, and in the Galileon model this can very 
easily lead to several numerical problems related with cases 
where the background evolution could develop 'fake' ghosts 
or Laplace instabilities |^ Consequently, although one could 
in principle work its way around these problems, for simplic- 
ity, we choose to sample h directly, even if this happens at the 
cost of having slightly slower runs. 

The trial-and-error search for C2 can be very time consum- 
ing and to speed it up we have adopted adaptive step lengths 
in setting trial values of C2- If the step length becomes very 
large, it is possible that the trial value of C2 becomes too far 
away from the right value so that the numerical solver of the 
EOM fails. However, we have checked that this is very rare 
with our algorithm, and by comparing with the results from 
fixing the step length we have confirmed that it does not affect 
our numerical constraints and conclusions. 



B. Scaling degeneracy 

As discussed in the previous section, one of the Galileon pa- 
rameters should be fixed when running the chains to break the 
scaling degeneracy of the Galileon equations. However, in or- 
der to gain some insight into how the scaling degeneracy man- 



^ 'Fake' here is in the sense that such cases could have been free of instabil- 
ities if C2 took the correct value for the consistent background evolution. 
However, this can only be done by simultaneously searching for C2 and h 
using trial-and-error and it is, therefore, much more time consuming. 
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FIG. 1. (Color online) Points accepted by the Metropolis-Hastings algorithm, including those sampled during the bum-in period, using 
WMAP9+SNLS-FBAO (blue dots). The points are projected onto different planes of the Galileon subset of the parameter space. All of the 
Galileon parameters were allowed to vary in order to manifest the scaling degeneracy. The remaining cosmological parameters were also 
sampled (not shown; see Fig.[3](. There are approximately 48000 points (around 8000 per chain) and the large red circles represent the fifty 
best-fitting points. 



ifests itself in the parameter space, we have first run a set of 
chains where all of the parameters are free to vary. The result 
is shown in Fig.[T| where we plot all the points accepted by the 
MCMC algorithm (cyan dots) for four different planes of the 
Galileon parameter subspace for chains constrained using the 
combined WMAP9h-SNLSh-BAO dataset. There are approxi- 
mately 48000 points (8000 per chain). We highlight the points 
with the fifty highest likelihood values (red circles). The rest 
of the cosmological parameters were allowed to vary as well 
(not shown in this plot, see Fig. [3] below). The parameters 
C3, C4 and C5 were sampled from the interval [—1000, 1000] 
to prevent the chains from spending too much time searching 
larger and larger values of the c„. The starting point of the 
Cn for each of the chains was set to be sufficiently close to 
zero so that the algorithm could quickly select the signs for 
the parameters which best fit the data. For practical reasons, 
we sample only positive (or only negative) values of (fi, since 
its absolute value can be very small. We point out, however, 
that this does not mean we are rulling out the regions of the 
parameter space where (fi < 0, since such regions can always 
be found by simultaneously fliping the signs of C3, C5 and (pi, 



as allowed by the scaling transformations of Eq. ( 22 1 



As expected from the scaling relations of Eqs. (22 1, one 
finds a long and naiTow region of degeneracy in the Galileon 
subset of the parameter space, along which the likelihood is 
kept constant. Note that although it seems that the best-fitting 



points are confined to C2 > —200, C3 < 400 and C4 > —600, 
this happens only because the points have reached the prior 
range limit of 1000 in the C5 direction, which therefore 'artifi- 
cially' constrains the other parameters. We have checked that 
the degeneracy region keeps increasing on increasing the size 
of the prior ranges. The difference in = — 21ogP (where P 
is the posterior) between the best-fitting and the fiftieth best- 
fitting points is Axit;, 50^/^ ~ — 1, but the likelihood does not 
change mono tonic ally along any direction of the parameter 
space. The result shown in Fig. [T|is in partial disagreement 
with the conclusions drawn in |72|. In the latter, the authors 
found a long and narrow region of degeneracy along which 
the likelihood decreases for values of c„ much larger or much 
smaller than unity. We agree that the long region of degener- 
acy exists. However, the likelihood does not change apprecia- 
bly towards larger values of the parameters c„, which is what 
one would expect in light of the scaling relation described by 
Eqs. ([22]). 



We see that the best-fitting points all lie in the region of 
parameter space where C2 < 0, C3 > 0, C4 < 0, C5 > 
when <fii > 0. This means that the sign of the fixed parameter 
C3 should be the same as (pi. For instance, if we were to fix a 
negative value of C3 while sampling only positive values of (pi, 
we would be discarding, a priori, the portion of the parameter 
space that contains the best-fittings points (red dots). Recall 
that the reflection symmetry only holds if the three parameters 



10 



C3, C5 and (pi all flip their signs. Throughout, the results will 
refer to sets of chains where C3 = 10 and ipi > 0, and we shall 
quote the final constraints in terms of the invariant combina- 



tions of Eq. (23 I. We stress that these combinations are also 



invariant under the simultaneous change of the sign of the C3, 
C5 and (pi parameters. In what follows we will refer to Cp and 

Cpl(?J'^ interchangeably, where p = 2, 4, 5. 



C. Parameter space constraints 

In Fig. [2] we show the marginalized two-dimensional like- 
lihood distributions and the corresponding 68% and 95% 
contour limits of the Galileon sector of the parameter 
space using WMAP9 data (top panels) and the combined 
WMAP9H-SNLSH-BAO data (bottom panels). Figure |3] shows 
the same but for the cosmological sector of the parame- 
ter space and the dashed black lines correspond to the 68% 
and 95% confidence limits for constraints on the ACDM 
model. We also show in Fig. |4] the marginalized one- 
dimensional distributions for constraints using the combined 
WMAP9h-SNLSh-BAO data. We point out that the constraints 
we derive for the ACDM model are in good agreement with 
those presented in the WMAP 9yr final result paper |7| (c.f. 
Tables|I]and|II|. 

As we discussed above, fixing one of the Galileon parame- 
ters breaks the scaling degeneracy of the Galileon sector of the 
parameter space, and therefore, we see in Fig. |2] that sensible 
constraints can now be derived. In particular, the parameters 



C2/C3 



2/3 



, C4/C3 



g^'^ and C5 jcj'^ are constrained to lie in a small 
region of the parameter space, very close to the theoretically 
unstable portion of the parameter space (see Fig. [TjbelowQ 
The shape of the contours shows also that these parameters 
are correlated to some extent. In Fig. [4j we see that the dis- 
tribution of Pip^ilpm,i is essentially flat, up to some statistical 
noise, but it decays very sharply for pip^i/pm^i > —6. This 
indicates that the upper bound of this parameter can play an 
important role in the physics of the Galileon |74|, as we dis- 
cuss below. 

The constraints on the Galileon sector of the parameter 
space do not change substantially when adding the SNLS and 
BAO data to the WMAP9 data. However, the same does not 
happen in the cosmological sector, in which case the preferred 
values differ quite significantly between the two data combi- 
nations, as shown in Fig. [3] In particular, when adding the 
SNIa and BAO constraints, the Galileon model prefers higher 
values of the total matter density ilmh^ ~ 0.145, as opposed 
to Qrnh^ - 0. 135 for WMAP9 alone (see Tables|l]and|ll{. The 
amplitude of the primordial fluctuations, Ag, also increases 
going from WMAP9 only to the combined dataset, but the 
preferred values of the expansion rate h, of the spectral index 
Us and of the optical depth to reionization, r, become smaller 



In I71II721 , the authors have also found that the high-likelihood regions 
tend to lie close to the regions associated with ghosts or Laplace instabili- 
ties. 



It is also interesting to compare how the constraints of the 
Galileon and ACDM models respond to different datasets. 
Here, we note that when using only the WMAP9 data, the 
Galileon and the ACDM models prefer more or less the same 
regions of parameter space with the exception of the Hubble 
rate h which is larger in the Galileon model, by just over la. 
However, adding the SNIa and the BAO data leads to devia- 
tions in all of the cosmological parameters whose marginal- 
ized two-dimensional distributions can differ by more than 
2(7. More specifically, when using the combined dataset, one 
sees in Fig. [3] that the Galileon model prefers more matter 
in the Universe, despite preferring slightly less baryons, and 
higher values of the parameters h and Ag relative to ACDM. 
The parameters Ug and t tend to be smaller in the Galileon 
model than in ACDM. The discrepancy between the con- 
straints of the Galileon and ACDM models for these cosmo- 
logical parameters highlights the importance of allowing all of 
the cosmological parameters to vary when extracting numeri- 
cal constraints, instead of fixing them to their best fit values in 
ACDM. For instance, note that the best-fitting values of h to 
the combined WMAP9h-SNLSh-BAO dataset are closer in the 
Galileon model than in ACDM to the independent determina- 
tions of the present-day expansion rate /i — 0.742 ±0.036 |8 j, 
h = 0.744 ± 0.025 t9J and ft = 0.743 ± 2.1 LIQJ. 

The Galileon model is able to fit the full WMAP9 data bet- 
ter than ACDM with a difference in of Ax^MApg ~ — 1-8 
(see Fig|6j negative values of Ax^ indicate the Galileon is fa- 
vored). On the other hand, once the SN and the BAO data 
are added, then the Galileon becomes less favored with a dif- 
ference of Ax^^MApg+SNLS+BAO ~ 8.6. This unveils a ten- 
sion between the CMB data and the geometrical tests in the 
Galileon model and it is worth trying to understand the rea- 
sons behind it. 

It is non-trivial that the constraints on the Galileon param- 
eters barely change between the two data combinations. In 
Il68ll69ll72l l74l. it was shown that the observations that only 
probe the background evolution of the Universe typically pre- 
fer larger values of the initial density of the Galileon field be- 
cause these minimize the strength and duration of the phantom 
evolution in the near past, therefore causing w to approach the 
de Sitter attractor w = — 1 sooner However, Fig|2] shows 
that adding the SNLS and the BAO data to the WMAP9 data 
does not lead to a significant increase in the upper bound of 
P<p.i/Pm,i (see Table |l]i. This is because the increase of the 
initial density has a larger impact on the time evolution of the 
gravitational potential than on the evolution of the background 
expansion rate. In |74| we showed that if the initial density of 
the field is too large, p,^,i/pm,i ^ 10^^ at Zi ~ 10^, then 
this generally leads to strong time variations of the gravita- 
tional potential and hence to an unacceptably large ISW effect. 
This is corroborated by the very sharp decrease in the poste- 
rior distribution towards large values of p^p^i/pm,i, as seen in 
Figs.[2]and[4] As a result, the ISW effect constrains the upper 
bound of this parameter so well that the constraint is barely 
changed when the background data is also considered. As we 
will see below, the parameters C4 and C5 are also very tightly 
constrained by the CMB data (c.f. Fig|7]|. 

In L74J . we have also shown that, if p^^ilpm.i is sufficiently 
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FIG. 2. (Color online) Marginalized two-dimensional posterior distributions and respective 68% and 95% contour limits obtained for the 
Galileon sector of the parameter space with WMAP9 data alone (top panels) and the combined WMAP9+SNLS+BAO datasets (bottom 
panels). The shading represents the distribution, where darker regions mean higher probability density. The posterior distribution and the 
respective contours were smoothed using a Gaussian filter that did not change the underlying results. In these chains the parameter C3 was held 
fixed at C3 = 10 and (pi was allowed to take only positive values. 



small, then different values will have no visible impact on the 
CMB power spectrum. The high-redshift evolution of w will 
still be sensitive to p^p^i/pm,i ^ lO"'^, but this dependency 
does not propagate to the expansion rate, which is only sensi- 
tive to w when the dark energy becomes the dominant com- 
ponent at low redshift. This is what makes the distribution of 
Pip,i/Pm,i, effectively, unbounded frombelowj^ as one can see 
in Figg]' 

Since the Galileon parameters are very well constrained by 
the ISW effect, the main consequence of adding the SNIa and 
BAO data to the CMB data would be to "tune" the total matter 
density in the Universe (which increases) and the present-day 
expansion rate (which decreases). This modifies the time evo- 
lution of H{a), but in a way that preserves the positions of 
the acoustic peaks of the CMB. This can be seen by com- 
paring the values for the age of the Universe in Tables [l] and 
[m which barely change between the two data combinations. 
However, these changes will have an impact on the ampli- 
tude of the CMB power spectrum, which causes the remain- 
ing cosmological parameters to explore the possible degen- 
eracies in order to keep as good a fit. As a result, they deviate 
from their best-fitting values obtained using the WMAP9 data 



^ Keeping in mind the positive sign of the energy density. Sufficiently nega- 
tive values of the energy density are typically associated with ghosts. 



alone. This "compensation" is, however, not perfect and the 
model is forced to reach a compromise between the regions 
of parameters favored by WMAP9 and SNLSh-BAO. Overall, 
the net result is a poorer fit to the combined dataset. This 
tension between the background constraints and the ISW ef- 
fect is similar to that highlighted in Refs. [69. ,71 J . However, 
here we want to emphasize that this is a tension between the 
cosmic background evolution and the cumulative effect of the 
time variation of the tensing potential, instead of one between 
the background evolution and the summed effect of the New- 
tonian potential^ In addition, this tension shows up at the 
largest length scales relevant for the ISW effect, where we are 
confident that the Vainsthein screening plays a negligible role 
so that there is no concern about the validity of linear pertur- 
bation theory. 



* In the conformal Newtonian gauge with the line element given by ds^ = 
[(1 + 2<I')dT2 - (1 - 2<I')dx2], the ISW effect probes ('i' + 'i>)/2 = 
<f), where 4> is the tensing (Weyl) potential, while particle dynamics, and 
therefore the growth rate, is sensitive to 'J'. Because of the nonzero 
anisotropic stress, Eq. (TTJ, these two potentials can be different in the 
Galileon model. 
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FIG. 3. (Color online) Same as Fig. |2]but for the cosmological sector of the parameter space. The corresponding contours obtained for 
the ACDM model are also shown for comparison (dashed contours). The scalar amplitude at recombination As refers to a pivot scale k = 
0.02Mpc"\ 




log[10^°AJ logioK,/p„.i] cjc^'' c,/c^'' 

FIG. 4. (Color online) Marginalized one-dimensional distributions obtained for the combined WMAP9+SNLS+BA0 dataset (solid blue). The 
distributions obtained for the ACDM model are also shown for comparison (dashed black curves). The scalar amplitude at recombination As 
refers to a pivot scale k — 0.02Mpc^^. In these chains the parameter C3 was held fixed at C3 — 10 and ipi was allowed to take only positive 
values. 
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FIG. 5. (Color online) Time evolution of the ratio of the Hubble ex- 
pansion rates of the Galileon and ACDM models, H / Hacdm (top 
panel) and of the Galileon field equation of state parameter w (bot- 
tom panel) for the maximum likelihood points of the chains using 
the WMAP9 data alone (dashed green) and the combined dataset 
WMAP9-l-SNLS-^BA0 (solid blue). The ACDM model used in the 
ratio of the upper panel is the best-fitting model obtained using the 
WMAP9-l-SNLS-fBA0 combined dataset. 



D. Cosmology of the best lit 

1. Background evolution 

In Fig.|5 we show the time evolution of the ratio H/H\cdm 
and w of the best-fitting points (shown in Table [llji obtained 
for the constraints with WMAP9 data alone (dashed green) 
and the combined WMAP9h-SNLSh-BA0 data (solid blue). 
Here, the ACDM model is the best-fitting model to the com- 
bined WMAP9h-SNLSh-BAO data. The expansion rate of the 
two best-fitting Galileon models agrees with ACDM during 
the radiation dominated era since these models all have the 
same energy density of radiation. However, in the matter era, 
the best-fitting Galileon model to the combined dataset has an 
expansion rate which can be up to w 4% higher than both the 
best fit Galileon model of the WMAP9 data and the ACDM 
model, as consequence of the also higher total matter den- 
sity (see Fig. [3] and Tables [11 and [ll|. At the end of the matter 
era, the expansion rates of the two Galileon models first de- 
crease to be approximately 4% smaller than ACDM, and then 
start increasing towards the present-day, at a « 0.5. This late 
growth is more pronounced in the former since the value of 
the present-day expansion rate, h, is larger 

In the best-fitting Galileon model to the combined 
WMAP9h-SNLSh-BAO dataset, the phantom evolution w < 
— 1 starts closer to the present-day and is also less pronounced 
than in the best-fitting Galileon model to the WMAP9 
data alone. This is mainly due to the larger value of 
logio [Pv'.i/Pm.i] in the former relative to the latter (c.f Ta- 
ble |n]i. As we showed in Fig. 1 of Q4\ (see also 1691), the 
lower Ptp^i/pm.i the more negative the values of w will be. 



This is because lower values of the Galileon density in the 
past will force the energy density to grow more drastically 
{w < —1) closer to today when the field starts to be driven 
towards the de Sitter attractor evolution. Recall, however, that 
the low redshift background data has little power to increase 
the initial Galileon density, since this parameter is very tightly 
constrained by the ISW effect (c.f. Fig.[2]i. 

It is interesting to compare the typical time evolution of the 
Galileon equation of state parameter with those reconstructed 
from the observational data in 190 1 (see also llllOlllll l). By 
using non-parametric Bayesian statistical techniques these au- 
thors have found that the currently available data favors dy- 
namical dark energy models. In particular, the preferred evo- 
lution of w is such that it crosses the 'phantom' line at z ^ 0.5 
from u; > — 1, and approaches w = — 1 from below at around 
the present time, which is not too different from what happens 
in the evolution of the background in the best-fitting Galileon 
model. The interested reader might wish to compare Fig.|5]of 
this paper with Fig. 1 of I.90J . 



2. Cosmic Microwave Background 

Fig.[6]shows the CMB temperature angular power spectrum 
of the best-fitting distribution values of the chains constrained 
using the WMAP9 data alone (solid green) and the combined 
WMAP9H-SNLSH-BAO data (solid blue) (c.f. Table The 
top right, bottom left and bottom right panels zoom into the 
three main regions of interests in the top left panel, which 
are, respectively, the low-l region which is most relevant to 
the ISW effect, the region of the first acoustic peak and the 
region of the higher peaks. For comparison, we also show 
the prediction of the best-fitting ACDM model constrained 
by the WMAP9 data alone, which we show as the open cir- 
cles with errorbars. We can see that the Galileon model gives 
a very good fit to the measured angular power of the CMB 
perturbations. For I > 100, the best-fitting Galileon mod- 
els give essentially the same prediction as ACDM, although 
the best-fitting Galileon model to the combined dataset (solid 
blue) predicts slightly more power in the interval of I that lies 
between the second and the third acoustic peaks. 

The most pronounced feature in the CMB predictions of 
the best-fitting Galileon models is, however, seen in the low- 
l region of the spectrum where the ISW effect dominates. 
In 1 74 1, we pointed out that in Galileon models the gravita- 
tional (Weyl) potential can undergo dramatic time variations, 
leading to a significant increase in the amplitude of the low-l 
CMB power. However, the results here show that the opposite 
is also possible, as it is the case of the best-fitting Galileon 
models. Indeed, even though the time variation of the grav- 
itational potential can be quite strong, it can also cancel out 
when integrated, yielding a weaker ISW effect for the best- 
fitting Galileon models, which even exhibit an inflection in 
the shape of the CMB power spectrum at Z « 6. This is in 
clear contrast with the ACDM curve, which keeps growing 
towards lower I, and such an interesting result can in principle 
be used to distinguish between the Galileon model and other 
dark energy or modified gravity models. It is well known that 
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FIG. 6. (Color online) CMB angular power spectra for the maximum likelihood points of the chains using the WMAP9 data alone (solid green) 
and the combined WMAP9+SNLS+BA0 dataset (solid blue). The power spectrum of the best-fitting ACDM model to the WMAP9 data alone 
is also shown for comparison (dashed black). The WMAP 9yr data is represented by the open circles with errorbars [7 |. The top right, bottom 
left and bottom right panels show zooms into the regions most relevant for the ISW effect, the first acoustic peak and the higher CMB peaks, 
respectively, of the top left panel. In the top left and top right panels at I — 4, from top to bottom, the curves correspond to ACDM, the best- 
fitting Galileon model using the combined WMAP-l-SNLS-l-BAO dataset and the best-fitting Galileon model using the WMAP9 data alone, 
respectively. In the bottom left panel, the ACDM curve is indistinguishable from the curve of the best-fitting Galileon model to the WMAP9 
data alone and at the CMB peak position, the curve of the best-fitting Galileon model to the combined WMAP9-I-SNLS-I-BA0 dataset is the 
one with lower amplitude. In the bottom right panel, the ACDM curve is indistinguishable from the curve of the best-fitting Galileon model to 
the WMAP9 data alone, but at the third CMB peak position (the second one shown in the panel), the curve with higher power is the curve of 
the best-fitting Galileon model to the combined dataset. In the rest of the bottom right panel the three curves are indistinguishable. 



the CMB signals at such low I are severely affected by cos- 
mic variance (note the large error bars on the WMAP data 
points), but nevertheless the CMB data does seem to prefer the 
Galileon models which exhibit such an inflection in the ISW 
power over those which do not. Note that this inflection is not 
a generic feature of the Galileon parameter space whose rich 
phenomenology allows it to have CMB power spectra which 
are also very similar to ACDM, as shown in 1741 . 

To see more clearly the effect of varying the c„ parame- 
ters, in the top right panel of Fig. [7] we plot the CMB power 
spectrum for the parameter points highlighted in the top left 



panel. In the latter, we zoom into a region of the C4 — C5 
plane constrained with the WMAP9h-SNLSh-BAO data (bot- 
tom right panel of Fig|2]l, to show the accepted points by the 
MCMC algorithm rather than the contours of the posterior dis- 
tribution. This enables us to resolve the fine details of the con- 
strained parameter space better than by evaluating the likeli- 
hood contour, which always requires some interpolation and 
smoothing. This is particularly useful in the Galileon model, 
whose parameter space can have a quite detailed shape. In- 
deed, we can see a gap between the accepted (cyan) and unsta- 
ble (grey) points, filled with low-likelihood points for which 



15 



the ISW effect is too strong. Furthermore, there is a 'hole' 
of low -likelihood points near the boundary of accepted points 
(where the point Bi is located), which is not visible in the 
contour representation. Points Ai, A2, A^, B and Bi are ob- 
tained from the best-fitting point A by keeping all the param- 
eters fixed except for C4 and C5, which are varied to explore 
the different parts of the C4 — C5 plane: Ai and A2 lie in the 
low-likelihood gap between the main sampled region and the 
unstable portion; Aj, is a typical accepted point; B is another 
accepted point which gives a fit nearly as good as A; and Bi is 
a rejected point in the low-likelihood 'hole' mentioned above. 

From Fig.|7]we can see that going from the best-fitting point 
A (green curve) inwards to the contour into another point still 
deep inside the 68% confidence limit (point A^, red line), we 
end up with nearly the same CMB power spectrum. On the 
other hand, moving outwards from the contour into points Ai 
(black line) and A2 (blue line) we get a sudden increase in 
the ISW power which is enough to render these models in- 
compatible with the data and, therefore, gives this region a 
low likelihood (in fact, not a single point was sampled by the 
MCMC search in the neighborhood of the points Ai and ^2). 
One sees that the ISW effect plays a key role in constraining 
the Galileon model, as already anticipated in our earlier work 
Il74l : the allowed region in the C4 — C5 plane corresponds to 
the values of these parameters for which the time evolution of 
the gravitational potential is mild enough that it yields small 
power on the largest angular scales (keeping in mind also that 
p^ i cannot be too large). The expansion rate evolution of the 
points in this region is essentially the same. The very strong 
sensitivity to the Galilean parameters can also be seen by look- 
ing at the CMB power spectra of points B (cyan line) and Bi 
(magenta line). Point B lies in a narrow layer of accepted 
points and gives an almost indistinguishable CMB spectrum 
compared to point A. Nevertheless, taking a small step to- 
wards the less favored region to point Bi leads to a significant 
increase in the power for I <11. Such an increase is however 
not as pronounced as that seen when going from point A to 
points Ai and A2, which is why some points in that region 
were still sampled during the MCMC search. 

The differences in the amplitude of the CMB spectra within 
the parameter space of the Galileon model show the advantage 
of using the full CMB data when constraining the Galileon 
model, over using only the WMAP distance priors as done in 
f68 , 7 1 72 1 . For instance, using only the information encoded 
in the positions of the acoustic peaks of the CMB would result 
in nearly the same fit for all the models shown in Figs.|6]and|7] 
although, in fact, these can be very easily distinguished and, 
in some cases, completely ruled out by the full CMB data. 



E. Future constraints 

The bottom left and bottom right panels of Fig. |7] show 
the power spectra for the linear matter density fluctuations 
and the weak lensing potential, respectively, corresponding 
to the models highlighted in the top left panel. In the bot- 
tom left we also show a recent estimate of the power spectrum 
of luminous red galaxies (LRGs) of the DR7 from SDSS UJ. 



The predictions of the Galileon power spectrum are shown 
at 2LRG = 0.31, which is the median redshift of the LRGs. 
This way one does not need to account for the growth factor 
to compare the theoretical spectrum with the observed one. 
However, there are at least three other effects that need to be 
considered, including baryonic bias effects, redshift space dis- 
tortions and nonlinearities. The latter includes the familiar 
mode coupling that develops between perturbations on differ- 
ent scales, but also, and this is particularly important in modi- 
fied gravity theories, possible screening effects (the Vainshtein 
screening for the case of Galileon models) which can render 
linear perturbation theory invalid on scales where it is usually 
taken to be a good approximation 1 251 l59l l97l l98l . We stress 
once again that this is the ultimate reason why we chose to 
leave growth rate and clustering measurements out of the con- 
straints in this paper, and conclude that a better understanding 
of the screening in Galileon models is necessary for reliable 
constraints using these datasets. 

However, in case the Vainshtein screening mechanism does 
not affect the linear predictions too much, a brief study of the 
linear matter power spectrum can still shed light on our under- 
standing of the model and indicate what we can expect when 
combining the different datasets. For this, one still needs 
to consider galaxy bias and redshift space distortion effects 
which link the theoretical matter power spectrum to the ob- 
served galaxy power spectrum. Galaxy bias arises because the 
baryons, compared with the dark matter particles, experience 
additional non-gravitational interactions which alter the spa- 
tial distribution of galaxies relative to dark matter Redshift 
space distortions result from the peculiar motions of galax- 
ies, which add to the Hubble flow, and result in systematic 
changes in the clustering power and pattern of the observed 
galaxies |96|. 

In the simplest cases, these effects are modeled by a 
scale-independent change in the amplitude of the spectra, 
though numerical simulations show that in reality these ef- 
fects are scale-dependent, especially for highly clustered ob- 
jects 112|. There is evidence which suggests that LRGs 
are biased tracers of the underlying dark matter distribution 
ll2l lll3tllT7l . which means that the observed galaxy spectrum 
should have a higher amplitude than the dark matter spec- 
trum predicted by linear theory, independently of the theory 
of gravity. This simple requirement already puts many of the 
Galileon models (e.g, most of those shown in Fig. |7]i at odds 
with the observations. The only exceptions are the points Ai 
and A2 whose linear matter power spectra are the closest to 
the ACDM prediction. In particular, point A2 predicts a lower 
clustering power than the standard ACDM model. However, 
these points are disfavored by the CMB data because they pro- 
duce too much ISW power (top right panel of Fig. [7]). This 
unveils yet another tension in the Galileon model, in addi- 
tion to that between the background evolution and the ISW 
effect and the one pointed out in ||72]| . In particular, one sees 
in Table |l] that the mean value of the variance of the power 
spectrum at fc = l/8ft.Mpc~^, cts, is erg « 1 for the Galileon 
model, as opposed to as ~ 0.8 for ACDM. With these re- 
sults we can anticipate that, providing linear theory holds on 
the scales of the observations, complementing our study with 
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FIG. 7. (Color online) The top left panel shows the points accepted by the MCMC algorithm (cyan dots) after the bum-in period for a zoomed 
region of the C4 — C5 plane constrained with the combined WMAP9+SNLS+BAO data. The highlighted points have the same parameters 
as the best-fitting Galileon model (point A, green circle) except for C4 and C5 which were chosen to explore other regions of the C4 — C5 
plane. The grey dots show the points sampled by the MCMC algorithm but that were rejected by being either associated with ghost or Laplace 
instabilities. The figure also shows the CMB power spectra (top right panel), linear matter power spectra at redshift ^lrg ~ 0.31 (bottom left 
panel) and angular power spectrum of the weak lensing potential (bottom right panel) for the Galileon models corresponding to the highlighted 
points in the top left panel and for the best-fitting ACDM model (dashed line) of the combined WMAP9-I-SNLS-I-BAO dataset. We also plot 
the WMAP 9yr data (open circles) |7| and the SDSS-DR7 LRG host halo power spectrum data (blue dots) |2| in the top right and bottom left 
panels, respectively, ^lrg is the mean redshift of the LRG sample. In the top right panel, at / = 2, for top to bottom the lines correspond, 
respectively, to the points A2 (goes over the plot limit), Ai, Bi, ACDM, A3, A and B. In the bottom left plot, from bottom to top, the lines 
correspond to point A2, ACDM, point Ai and are indistinguishable from there on. In the bottom left panel, from top to bottom, at the peak 
position of the ACDM curve (dashed black) the lines correspond to the points B, A, A3, Bi, Ai, A2 and ACDM, respectively. 



clustering measurements will result in even tighter constraints 
than those quoted in Tables [l| and [II] The only salvation would 
arise if any modifications introduced by the Vainshtein mech- 
anism are able to bring the clustering predictions closer to the 
observations. 

In 1741 we saw that the power spectrum of the weak lensing 
(Weyl) potential (p can also lead to very distinctive observa- 
tional signatures of the Galileon model. The weak lensing 
signal is the integrated effect of the Weyl gravitational poten- 
tial along the line-of-sight from today to the epoch of recom- 



bination 111 181 II 19l . In Fig. [7] we see that the angular power 
of the lensing potential in the Galileon model is higher than 
ACDM on all scales for all of the models shown, which in- 
cludes models whose CMB power spectrum fit the WMAP9 
data sufficiently well. This once again provides another po- 
tentially powerful way of discriminating between the Galileon 
and the standard ACDM model. Moreover, in the case of 
the weak lensing power spectrum, just as in the CMB case, 
the length scales are large enough for one to not worry about 
whether or not the Vainshtein effect has an impact on the va- 
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lidity of linear perturbation theory. However, the currently 
available data still has little constraining power, even for the 
ACDM model, and therefore no decisive conclusions can be 
drawn for now. Nevertheless, this shows that weak lensing 
measurements, if sufficiently accurate, may provide a promis- 
ing way of constraining the Galileon model, which can be in 
many aspects representative of other variations of the general 
Horndeski model. For this, it would be interesting to forecast 
how accurate future lensing measurements would need to be, 
like the ones of Euclid and Planck, in order to distinguish the 
GaUleon model from ACDM. 



V. DISCUSSIONS AND CONCLUSIONS 

We have studied and constrained, for the first time, 
the full parameter space of Galileon gravity models, 
which is composed of the six cosmological parameters 
ilco,^bo, h,T, Us, \og[lO^^ As]; plus the Galileon parame- 
ters, using the latest observational data on SNla, BAO and 
the full temperature anisotropy power spectrum of the CMB. 
This is the first time the full CMB data, rather than just the in- 
formation encoded in the positions of the acoustic peaks, has 
been used to place constraints on the Galileon model. Our re- 
sults were obtained using a modified version of the CAMB code, 
which solves for the cosmological background and linear per- 
turbation evolutions fT4\. We implemented these solutions in 
the publicly available CosmoMC code to sample over the pa- 
rameter space. 

We found that there is a set of scaling relations which can 
be used to change the Galileon parameters without changing 
the physics. We showed that when all of the parameters are 
allowed to vary, the scaling relations conspire to produce an 
unlimited degeneracy region in the Galileon parameter space 
along which the likelihood barely changes. A particularity of 
this degeneracy region is that it is symmetric in the signs of 
the parameters C3, C5 and (pi. A similar scaling relation was 
also found in f67l (and more recently in [71 1). However, in 
this paper we treat it in a way which allows us to directly con- 
strain the physical parameters of the Lagrangian. Due to the 
scaling degeneracy, the dimensionality of the Galileon param- 
eter space is reduced by one and we need only to constrain 
specific combinations of the parameters which are invariant 
under the scaling transformation. This result does not fully 
agree with that of 1721 . where the infinite degeneracy does not 
seem to manifest itself, despite, apparently, allowing all of the 
parameters to vary. The latter result is not compatible with the 
existence of the scaling degeneracy of the Gahleon equations. 

We tried various combinations of data to assess the ability 
of the Galileon model to fit different observations. In partic- 
ular, the Galileon model can fit the WMAP9 data better than 
ACDM with a difference in between the best-fitting models 
of Ax^ « — 1.8 (negative values indicate the Galileon model 
is favoured). However, when the CMB data is combined with 
measurements from SNla and BAO, the Galileon model be- 
comes less favored with Ax^ ~ 8.6. 

We found that the CMB data constrains the Galileon sub- 
space of parameters extremely tightly through the ISW ef- 



fect; this prediction depends sensitively on the exact values 
of the parameters. The best-fitting Galileon models are those 
for which the CMB power spectrum, in the ISW region, de- 
cays at low values of I. This result is in contrast to ACDM, 
where the ISW power increases at low I (c.f. Fig. |6]). We 
saw that the c„ parameters are mildly correlated with one an- 
other (c.f. Fig. [2]) and that tiny deviations from the best-fitting 
regions can quickly lead to a dramatic increase in the ISW 
power (c.f. Fig. [7]i, rendering the models incompatible with 
the data. Moreover, the ISW effect is also responsible for the 
very sharp cut-off in the posterior distribution of the initial 
energy density of the Galileon field above p^p.il Pm,i ^ 10^^ 
(c.f. Fig. [4]i. This dependency of the ISW effect on the 
Galileon parameters is so strong that adding the low redshift 
data of SNla and BAO barely changes the allowed regions of 
the parameter space (c.f. Fig.|2]). 

The constraints on the cosmological parameters in the 
Galileon model, on the other hand, show a much stronger de- 
pendence on the combinations of data used. Going from the 
constraints using WMAP9 data alone to the constraints using 
the combined WMAP9h-SNLSh-BAO dataset, we find that the 
best-fitting values of the total matter density, il,,„o/i^, and am- 
plitude of the primordial scalar fluctuations, Ag, shift towards 
larger values; while the expansion rate h, the scalar spectral 
index of the primordial power spectrum of scalar perturba- 
tions, Ug, and the optical depth to reioinization, r, become 
smaller The fact that the different datasets prefer different 
regions of the cosmological parameter space in the Galileon 
model shows the existence of a tension between the cosmic 
background evolution and the CMB data, which, ultimately, 
leads to a poorer fit to the combined datasets. 

We compared the constraints on the cosmological parame- 
ters in the Galileon model and in ACDM. These two models 
prefer, approximately, the same best-fitting parameters when 
using the WMAP9 data alone, with the exception being the ex- 
pansion rate h, which is higher in the Galileon case. However, 
in the case of the combined WMAP9h-SNLSh-BAO dataset, 
the best-fitting parameters of the two models are consider- 
ably different, with some of the marginalized two-dimensional 
posterior distributions being discrepant by more than 2a (see 
Fig. [3]l. In particular, the Galileon model predicts larger val- 
ues for fimoh^, h and As, but lower values for 71 <; and r, than 
ACDM. This shows the importance of varying all the cosmo- 
logical parameters when comparing to observations, since in 
this way a broader range of degeneracies in the parameters can 
be explored, leading therefore to more robust and informative 
conclusions. 

To try to anticipate how the Galileon parameter space can 
be further constrained by clustering data, we looked at the lin- 
ear matter power spectra of the best-fitting Galileon models 
to the combined WMAP9h-SNLSh-BAO dataset. The situa- 
tion here is more complicated given the uncertainties about 
the impact of the highly nonlinear Vainstein screening on the 
scales for which growth and clustering data are available. 
In particular, the major issue is whether or not linear the- 
ory is still valid on the relevant scales. Nevertheless, pro- 
vided that linear theory is indeed a good approximation, we 
have found that the best-fitting models using the combined 
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TABLE I. Mean values and uncertainties (la level) of the marginalized one-dimensional distributions of the full cosmological parameter space 
for the Galileon (except log^o [p^,i/pm,iQ and ACDM models using WMAP9 alone and the combined WMAP9+SNLS+BA0 data. The 
scalar amplitude at recombination refers to a pivot scale k — 0.02Mpc^^. The subscript "i" refers to quantities evaluated at z = Zi = 10^. 
We also show the mean values of the age of the Universe and of erg at redshift 2 = 0, which are derived parameters. 



Parameter Galileon Galileon ACDM ACDM 

(WMAP9) (WMAP9+SNLS+BA0) (WMAP9) (WMAP9+SNLS+BA0) 



SIM)''- 

h 



logio [10^° As] 

/ 2/3 
C2/C3' 

4/3 



C4/c: 

/ 5/3 
C5/C3' 



0.0227: 



r, ,,^+0.007 



0.817 ±0.04 
0.971 ±0.018 



0.0898t°:«°«^ 



o ,Q4+0.034 
J.iU4_0.033 



-5.18" 



4 21+0.39 

^■^^-0.38 



-0 ^(^\+°■^* 

U.iDi_o 3g 

042+°"^'' 



0.02181 ±0.00041 
0.1255 ± 0.0026 
0.735 ±0.012 
0.945 ±0.011 



0.0775 



+0.0052 
'0.0057 



q 1^7+0.025 
O.iy:(_o,o24 



-5.44" 



*-U*-0.34 



r, 1 71 +0.035 
J^-0.032 



046+" °" 



02259+"-'""''^'' 

r, 1 1 9+0.053 
'J-J^J^^-0.055 

0.705 ±0.027 
0.967 ±0.015 
0.0889l«:r6^ 

q 107+0 031 
O.IU ( _o.029 



0.02244 ± 0.00045 
0.1139 ± 0.0028 
0.695 ±0.012 
0.963 ±0.011 



0.0867l«:r6? 



3.112 ±0.028 



(78 (2 = 0) 0.938 ±0.038 0.988 ± 0.022 0.810 ± 0.026 0.817 ±0.020 

Age(Gyr) 13.49lo:l8 13.770to.o87 13.75 ±0.13 13.793 ±0.093 

Since the distribution of logjQ [p<^,i/ Pm.i] is unbounded from below, it does not make sense to quote its mean. Instead, we quote the highest value sampled 
by the chains to give an estimate where the distribution cuts-off. 



WMAP9+SNLS+BAO dataset predict too much clustering of 
matter on small scales, if 6 > 1 (where h is the bias factor). We 
did manage to find points in the parameter space which predict 
a linear matter power spectrum closer to ACDM, but the ISW 
effect for these models is, however, too strong to be compat- 
ible with WMAP. Moreover, the best-fitting Galileon models 
prefer values of erg which are substantially larger, erg « 1, than 
in ACDM, where erg ~ 0.8 (c.f. Table This suggests that 
there is a tension between the CMB and matter clustering for 
this model, which might be difficult to avoid so that adding 
clustering measurements may significantly tighten the con- 
straints. But again, because of the concerns discussed above, 
in order to have a cleaner and robust constraint, we prefer to 
leave the inclusion of clustering and growth data until we have 
a better understanding of the nonlinear effects in the Galileon 
model. 

We have also looked at the power spectrum of the lens- 
ing potential for the best-fitting models constrained by the 
WMAP9+SNLS+BAO dataset and found them to have, sys- 
tematically, more power than the ACDM prediction on all an- 
gular scales. The currently available weak lensing data is still 



not good enough to derive meaningful constraints and, there- 
fore, one cannot yet tell if another tension with data will arise 
from here. Nevertheless, the weak lensing signal seems to be 
able to provide a complementary way of distinguishing be- 
tween the Galileon and the ACDM model. Moreover, the an- 
gular scales relevant for the lensing potential power spectrum 
are large enough for one to assume that the Vainshtein screen- 
ing effects should be negligible. As a result, it would be in- 
teresting to forecast the necessary accuracy of future lensing 
observations to distinguish between the Galileon model and 
ACDM. 

In summary, we find that the currently available CMB and 
background data have the power to place very tight constraints 
on the Galileon parameter space. The model is slightly less 
favored than ACDM but one cannot yet safely rule it out us- 
ing data from the CMB and low redshift geometrical tests. 
On the other hand, the best-fitting Galileon models make very 
distinctive predictions for the values of fimo, h^Us, As,t and 
(Tg, which opens for the possibility of the Galileon model to be 
distinguished from ACDM by higher-significance cosmologi- 
cal data, or by currently available measurements of the galaxy 
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TABLE II. Maximum likelihood points in the MCMC chains together with the best-fitting value of ~ — 21ogP for the Galileon and ACDM 
models using WMAP9 alone and the combined WMAP9+SNLS+BA0 data. The scalar amplitude at recombination refers to a pivot scale 
k = 0.02Mpc~^. We also show the initial value of the Galileon field time derivative and of the age of the Universe. The subscript 'V' refers 
to quantities evaluated at z = Zi = 10^. 



Parameter 


Galileon 
(WMAP9) 


Galileon 
(WMAP9+SNLS+BAO) 


ACDM 
(WMAP9) 


ACDM 
(WMAP9+SNLS+BAO) 


x' 


7556.87 


7989.97 


7558.64 


7981.42 




0.02233 


0.02178 


0.02238 


0.02237 




0.116 


0.125 


0.113 


0.113 


h 


0.789 


0.735 


0.697 


0.693 




0.957 


0.947 


0.966 


0.963 


T 


0.0777 


0.0680 


0.0879 


0.0892 


login [10^°AJ 


3.112 


3.127 


3.117 


3.116 


loglO [Pv,i/Pm,i] 


-9.15 


-6.51 






1 2/3 
C2/C3' 


"3.67 


-3.59 






/ 4/3 
C4/C3' 


-0.195 


0.199 






/ 5/3 
C5/C3' 


0.0485 


0.0501 






- 1/3 

Age (Gyr) 


7.42 X lO"^'' 
13.572 


2.31 X 10"" 
13.778 


13.789 


13.819 



clustering, provided linear theory is applicable. 
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